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Abstract 

We review multifragmentation data obtained at the SIS/GSI accelerator using heavy 
ion beams with (O.l-l)A GeV together with the ALADIN and FOPI experimental 
setups. 



1 Introduction 

Multifragmentation, that is the emission of several intermediate mass frag- 
ments from a hot nucleus, is a phenomenon observed in nuclear reactions, 
using light and heavy projectiles over a wide incident enery range. Besides the 
scientific curiosity to understand what happens, some of the hopes in these 
studies are that one can learn more about the tendency of fermionic nuclear 
matter to appear in clusters and, perhaps eventually, about the topology of the 
nuclear phase diagram, in particular the evasive liquid-to-gas transition and, 
from the explosive features of some of the reactions, about nuclear matter com- 
pressibility. Such information is of high interest in other fields as well, such as 
astrophysics. As the underlying mechanisms are likely to be rather complex, 
statistical approaches have been invoked from the beginning in order to try to 
understand the experimental data. Evidence for non-equilibrium effects have 
lead to the parallel development of sophisticated transport models as well. 

As this written contribution represents essentially the contents of a one-hour 
lecture, it is highly limited in scope: I describe experimental multifragmenta- 
tion data measured in the past years at GSI, Darmstadt (Germany). The many 
facets of multifragmentation are studied in laboratories all over the world, how- 
ever. I refer the reader to a recent workshop on Multifragmentation for more 
encompassing information on this subject. 

I shall start by describing the ALADIN and FOPI experimental setups and 
the way observed multiparticle events are sorted. Then I review fragment yield 
distributions in both central and peripheral collisions together with attempts 
to derive from them apparent chemical temperatures. Detailed information on 
momentum space distributions will allow me to discuss subjects such as the 
determination of the caloric curve, the influence of flow and the system-size 
dependence. Finally, an isospin tracer experiment, performed to study the 
degree of equilibration, will be described. 

^ E-mail address: W.Reisdorf@gsi.de 
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2 Apparatus 



Two large detection systems have been set up at the heavy ion synchrotron 
SIS at GSI Darmstadt, ALADIN (Fig. |I|) and FOPI (Fig. g). 
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Figure 1: Schematic view of the ALADIN spectrometer in the bending plane 
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Figure 2: Schematic view of the FOPI detector. 



The ALADIN spectrometer was designed primarily to look at fast 
projectile-like spectator fragments focussed into a forward cone of about ib5°. 
The core of the apparatus consists of a magnet operated at a bending power of 
1.4 Tm, a time projection multiple-sampling ionization chamber TP-MUSIC 
and a time-of-flight (TOF) wall. In some experiments complementary detectors 
are installed, such as a large area neutron detector LAND that identifies neu- 
trons emitted close to 0°, and a series of multidetector hodoscopes consisting 
of Si-Cs(Tl) telescopes, CaF2 plastic phoswich detectors and Si-strip detectors 
arranged so as to cover angles between 6 — 58°. One of the purposes of the 
TP-MUSIC detector is to allow precise tracking of charged particles passing 
through its active volume. The version MUSIC-II P] is shown in Fig. |[ 




time 



Figure 3: Illustration of the operation of 
the MUSIC II detector From ref. § 

It consists of three active volumes with the electric drift fields in adjacent 
sections perpendicular to each other, two for the measurement of the horizontal 
and one for that of the vertical position and angle of the particle track. To 
allow multiple sampling of the particle signals, each anode is subdivided into 16 
stripes with a width of 3 cm each. The anode signals are recorded using flash 
analog-to-digital converters (ADCs) with a sampling rate of 16 MHz. Better 
than unit resolution for nuclear charges from Z=8 all the way to Au (Z=79) 
is obtained (Fig. Charges between Z=2 and 15 can be resolved with use of 
the scintillator strips in the TOF. 

The FOPI apparatus was built to study central heavy ion reactions in the 
energy range 0.1-2 A GeV. Particle identification, again, is done by a com- 
bined energy-loss, time-of-flight and magnetic rigidity analysis. The magnet 
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is a superconducting solenoid operated at 0.6 T (see Fig. 0). Time of flight 
and energy-loss are determined by about 1000 scintillator detectors arranged 
octogonally in the downstream part of the detector (PLASTIC WALL) and 
as a barrel inside the magnet. At lower incident energies (E/A < 400 MeV) 
a set of gas ionization chambers is inserted in front of the PLASTIC WALL 
to allow identiflcation of slower heavy clusters (up to Z about 15). Inside the 
magnet, the particles are tracked in drift chambers, the Central Drift Cham- 
ber, CDC, and the HELITRON (usually installed in front of the CDC, but not 
shown here). Fig. ^ allows to visualize individual particle tracks in a speciflc 
event determined with use of local or global (Hough transform) track flnding 
methods. 




Figure 5: Event display in the Central Drift Chamber. A cross section in the plane 
perpendicular to the beam axis is shown 

By matching these tracks to the outer barrel scintillators one obtains the 
energy loss (from the ionization in the CDC gas), the track curvature and the 
time-of-flight (from the Barrel ). Fig. ^ shows well separated pions , protons, 
deuterons and tritons. 
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Figure 6: Identification of charge-one par- 
ticles in the velocity versus magnetic rigid- 
ity plane 



Figure 7: Evolution of a heavy ion collision 
spectators and participants 



3 Event sorting 

Heavy ion collisions lead to events with 
very different topologies. It is highly 
desirable to sort these events according 
to some criteria, such as the impact pa- 
rameter h. As h itself cannot be mea- 
sured directly one tries to find an ob- 
servable that is believed to be strongly 
correlated to h. The methods used by 
both ALADIN and FOPI are easiest 
to visualize in the participant-spectator 
picture illustrated in Fig, ^. In experi- 
ments with the ALADIN spectrometer one selects impact parameters by mea- 
suring the size of the projectile spectator: one adds up the charges of all the 
fragments with Z > 1 moving forwards in a narrow rapidity window, see Fig. p. 
This sum has been dubbed Zhound- large Zbound obviously correspond to rather 
peripheral collisions. In experiments with the FOPI detector which are pri- 
marily performed to study fireball (participant) physics, a standard method 
is to determine the charged-particle multiplicity. The underlying idea is that 
'fireball matter' emits more particles than 'spectator matter'. In the early 
phase of running the detector the multiplicity, PM, observed in the PLASTIC 
WALL, which covers only laboratory angles forwards of 30°, was taken. Due 
to this geometric limitation, overcome later when the full detector became 
available, other more powerful methods to 'maximize' the fireball (i.e. choos- 
ing very central collisions) were developped. Fig. § illustrates the increasingly 
'compact' momentum space topology of Z=4 fragments emitted in Au on Au 
reactions at 2 50 A MeV as selection is done, first, according to high Plastic Wall 



5 



multiplicities PM (lower panel, where still spectator peaks are visible), second, 
according to large 'ERAT', the ratio of total transverse to longitudinal kinetic 
energies (center panel) and finally, using large ERAT and low directivity, i.e. 
high azimuthal symmetry. The directivity is obtained from Yl'^Zut/ Z\ut\ 
{uj = ±1 forward/backward, Z fragment charge, Ut transverse 4- velocity). In 
the figure both the rapidity (abscissa) and the transverse 4-velocity (ordinate) 
are scaled with the respective projectile values, placing the midrapidity value 
at zero. The ERAT and PM selections correspond to cross sections of 200 mb, 
the additional low-directivity condition cuts the cross section down to 70 mb. 
For more details see Pl. 
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Figure 8: (a): Rapidity spectra measured 
in the reaction Au+Au at lOOOA MeV for 
fragments with Z=2,3,4 and 6. (b): Ra- 
pidity spectra for Z=2 fragments in central 
collisions and (c): in peripheral collisions. 
From ref. H 



Figure 9: Invariant cross section contour 
plots (scaled units) for Be (Z=4) fragments 
observed in the system Au+Au at 250A 
MeV. The three panels correspond to differ- 
ent event selections, see text. From rcf. H 
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4 Cluster production: peripheral and central collisions. 
Universality. 

Multifragment decays of projectile-spectator sources have been studied exten- 
sively with the ALADIN spectrometer. [@] - One of the hopes was that spec- 
tator matter would be 'gently' excited, initially uncompressed nuclear matter, 
that would equilibrate relatively fast, then expand slowly under heat pressure 
and enter into the unstable spinodal regime where its decay would show fea- 
tures of a liquid-to-vapour transition. Evidence that multifragmentation was 
a phenomenon following complete equilibration was obtained by the observed 
'universality' of the decay pattern. As shown in Fig. |T0|, if one plots the ob- 
served mean multiplicity of intermediate mass fragments (IMF, Z = 3 — 30) 
against the size of the spectator charge , characterized by Zbound, (i-e. against 
the impact parameter for a fixed target- projectile system), one finds a depen- 
dence that remains invariant against the variation of the incident energy. 0] 




Figure 10: (a): Mean multiplicity of in- 
termediate mass fragments as a func- 
tion of Zjyound for the reaction Au+Au 
at E/A = 400, 600, 800 and 1000 
MeV. g 



Figure 11: Left panel: Mean multiplicity of 
intermediate mass fragments as a function of 
Zbound for the reactions ^^^U on ^^"^Au (circles), 
^^^Au on ^^^Au (squares) and ^^^Xe on ^^^Au 
(triangles) at E/A= 600 MeV. Right panel: 
the same data after normalizing both quanti- 
ties with respect to the atomic number Zj, of 
the projectile, g] 



Other observables characterizing the population of partition space were 
found to be 'universal' as well. Provided the projectile beam (Au) was not 
changed, universality was not lost when the Au target was replaced by lighter 
target materials, all the way down to Be. Further 'ALADIN plots', obtained 
by fixing the target to Au, but varying this time the size of the projectile 
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from U down to Xe (left panel Fig. |Tl]), are again merged into one curve when 
scaling the IMF multiplicity with the projectile size (right panel Fig. |TT]). 

The independence of fragment charge distributions on the size of the source 
could also be observed in central fusion-like reactions at lower incident energies 
using the Air multidetector INDRA at the GANIL accelerator, Caen, France, 



see Fig. |g. [g Studying the systems ^^sXe+Sn at 32A MeV and ^^^U+^^^Gd 
at 36A MeV, it was found that the charge distributions for fragments with 
charge Z > 4 were essentially identical, if scaled with the system size, or, 
as done in the figure, with the average multiplicity of charged particles with 
Z > 4 (the ratio of the latter is close to the ratio of the system sizes). It was 
argued in ref. |10] that the available excitation energy, 7A MeV, was the same 



for both systems. 
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Figure 12: Charge distributions in the reactions Gd+U at 36 A MeV (circles) and Xe+Sn 
at 32A MeV (triangles). Insert: Multiplicity distributions for fragments with Z > 5 in the 
same systems. After ref. 11^] 



Beaulieu et al. [O pointed out earlier that the excitation energy per nu- 
cleon was the relevant parameter for sorting out multifragmentation data. In 
Fig. |13] the authors show evidence for scaling of system-size normalized IMF 
multiplicity data, when plotting them against the excitation energy per nu- 
cleon: quasi-projectile decay data for the reactions ^^Cl+^^''Au at 43A MeV 
and ™Ge+Ti at 35A MeV are shown together with ALADIN spectator decay 
data in the reaction Au-|-Au at 600A MeV. The non-trivial determination of 
the excitation energy deposited into the presumably thermalized part of the 
system will be discussed in the next section. As pointed out in ref. ||Tl]| there 
is an important difference between 'quasi-projectile' data (E/A < 50 MeV) 
and 'projectile-spectatator' data (E/A > 600 MeV): in the former case the 
average mass of the emitting sources is almost constant due to total-charge 
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Figure 13: Average number of intermedi- 
ate mass fragments (scaled by the system 
size) as a function of the excitation energy 
per nucleon. Triangles: "^^Cl + ^^^Au at 
43A MeV, squares: ™Ge-l-Ti at 35A MeV, 
crosses: ALADIN data. After ref. ||ll| 
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requirements in the experiment, whereas in the latter case the Au projectile 
spectator decreases from mass A ~ 190 to A ~ 50 over the excitation energy 
range in the figure. 

Before closing this section, let us 
briefly present the 'fireball' part of the 
story. In very central Au on Au col- 
lisions above lOOA MeV one finds al- 
most perfectly exponential charge dis- 
tributions (Fig. 0), the slope merely 
becoming steeper as the excitation en- 
ergy is raised. As shown in the Fig. |1^, 
the probability for a proton to appear 
attached to at least one other nucleon 
decreases rather slowly from about 80% 
at 200A MeV to still 50% at lA GeV. 
As in these energies the 'available' 
c.o.m. energies are well beyond typi- 
cal nucleon binding energies, dynamic 
mechanisms, such as cooling by adi- 
abatic expansion, or non-equilibrium 
processes must be invoked to ex- 
plain such high degrees of clusteriza- 
tion. 




2 4 6 8 10 
nuclear charge Z 



Figure 14: Charge distributions in central 
collisions of Au on Au at 150, 250 and 400A 
MeV. il 



5 Chemical temperatures 

Frequently 'Statistical Multifragmentation Models' (SMM) are used to com- 
pare with data under the assumption of at least chemical equilibrium. The 
physics input involved are degrees of freedom of hot liquid drops modeled in 
analogy to the Weizsacker formula for cold nuclei, a more or less extended list 
of known discrete levels and some global level density formula to extrapolate 
to unresolved continuum excited states. Subtle sampling methods (within the 
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Figure 15: Degree of clusterization (%) of protons in central Au on Au collisions. 



grand-canonical, canonical or even microcanonical framework) allow one then 
to dial typical 'events' under the minimal assumption of basically homoge- 
neous available phase space population. Before being able to compare with 
experiment the 'primordial' fragment distribution has to be corrected for later 
sequential decays and the adiabatic pickup of potential field energy prevalent 
at the moment of 'freeze out'. (Usually only the Coulomb part is taken into 
account). Comparison of one such model [jl3|, Q with the ALADIN data is 
shown in Figs, [l^, |l^, |18| taken out of the thesis work ||^ of T. Odeh. Typically, 
one assumes freeze out at a nucleon density p = l/3po (Po is the saturation 
density), and adjusts the excitation energy (which such static theories cannot 
predict) until the data are reproduced. As the ALADIN curve must come 
down to zero at the left end because of reaching zero spectator size, as well as 
at the right end (high values of Zbound) because grazing collisions do not lead to 
any sizeable excitation, there is sensitivity to the theory, and in particular to 
the assumed excitation energy, only for intermediate Z^ound values. Excitation 
energies of about 4, resp. 10 MeV per nucleon are deduced for spectator nuclei 
with masses 160, resp. 40. Once the ALADIN curve is reproduced, (upper 
panel Fig. ^), one finds that other observables characterizing partition space 
are reproduced as well: see for instance the lower panel of Fig. ^ where the 
asymmetry parameter ai2 = (^max - ^2)/(^maa; + ^2) is plotted {Zmax and Z2 
are the largest and next-to- largest fragment charges). 

The chemical temperature deduced from the model is displayed in Fig. ^ 
and turns out to be rather flat in the Zbound range where there is sensitiv- 
ity to the data adjustment. A well known method |15| to obtain density- 
independent temperatures consists in evaluating double isotope yield ratios, 
such as (^Li/^Li) / ('^He/^He). The deduced temperature, dubbed T(HeLi), 
corrected for sequential decays, is also plotted in Fig. |l^ together with that 
deduced from the statistical model applying to it the same treatment as to the 
isotope data. The remarkable consistency of these apparent temperatures has 
to be confronted, however, with the failure of the model to predict correctly 
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Figure 16: Average intermediate mass frag- 
ment multiplicity (upper panel) and aver- 
age charge asymmetry (lower panel) as a 
function of Zbound for experiment (full dots) 
and statistical multfragmentation theory 
(open dots). The dotted and dashed curves 
show the model sensitivity to A ±15% vari- 
ation of the excitation energy. M 



Figure 17: The temperature T(HeLi) as a 
funtion of Zt,ound- Both projectile (600 and 
lOOOA MeV) and target (lOOOA MeV) spec- 
tator data are given. The dashed line rep- 
resents T(HeLi) as determined with use of 
the statistical model and the full line gives 
the model's breakup temperature. |p| 



the ratios of the hydrogen isotopes, Fig. 18 
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Figure 18: Proton to deuteron (left) and deuteron to triton ratios as function of Z\yound in 
the reaction Au+Au at lA GeV. The experimental data are represented by data, statistical 
model calculations are given by the smooth lines. From ref. H 
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Can the clusterization in the fireball, Figs. |M|, |T5|, also be understood in 
terms of chemical equilibrium? In view of the simple exponential shapes of 
the charge distributions, it does not come as a surprise that, given a fixed 
freeze-out density p, and the freedom to adjust the overall normalization, one 
can find an apparent chemical temperature T that reproduces the measured 
slopes well. The (p, T) pairs compatible with the FOPI data of Fig. in the 
framework of the Quantum Statistical Model |TB[, QSM, are plotted in Fig. 
adapted from ref. 
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Figure 19: Temperature versus freeze-out Figure 20: Temperature T(HeLi) (full sym- 
density as determined from the compar- ^o^^) excited state temperatures (open 
ison of the experimental charge distribu- symbols) as a function of incident aiergy 
tions with those calculated with QSM for P^r nucleon in Au on Au collisions. 
Au+Au at 150A, 250A and 400A MeV. 
Some isentropic curves (constant S/A) are 
also shown. ||l^ The shaded area represents 
the spinodal region. ||l^ Cross: inferred 
from adjacent Figure. See text. 

As the (p, T) pairs run approximately along lines of constant specific entropy 
(S/A), one can use this kind of study to deduce the latter. |]12| Just for the 



sake of orientation, the spinodal region, calculated for isospin symmetric 
matter, is also indicated in the figure. Unless one allows for rather high freeze- 
out densities (such as 0.8po), these apparent temperatures are rather low as 
we will see later. 

Fireball matter was also studied in a Au on Au experiment with the AL- 
ADIN setup. [0 The temperature T(HeLi) is found to rise steadily with in- 
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cident energy, Fig. and is compatible with the FOPI results in Fig. |T9| if 
the overlapping T(HeLi) point at 150A MeV (cross in Fig. ^) is inserted at a 
density close to 0.5po • 

Again, a warning is given by Nature: Temperatures deduced []^ from 
excited-state populations (^Li, "^He, ^Be) seem to cluster around 5 MeV and 
are invariant against variation of the incident energy. See Fig. |2D. 



6 Peripheral collisions: Kinetic energies, calorimetry 

The considerations in the previous section left out the information on momen- 
tum space distributions. 
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Figure 21: Energy spectra of light charged particles in three different bins of Zbound- Solid 
dots: Experimental data obtained at 150° in the laboratory. Open symbols: SMM calcula- 
tions, n 
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Figure |21| (light charged particles) and ^ (IMF's) summarize, for Au on 
Au at lA GeV, the remarkable differences between the experimentally ob- 
served kinetic energies and the predictions of the statistical approach, that 
was adjusted to render correctly the yields of Z > 2 spectator fragments. 
Taking a closer look at Fig. ^ (the normalization of theory was done at 
the ^He level, lowest panels), one notices that the deviations between the- 
ory and experiment increase from more peripheral (large Zi,ound, right col- 
umn) to more central collisions (left column) and decrease as the fragment 
mass is raised from A = 1 (top row) to A = A (bottom row). The av- 
erage kinetic energies of the IMF's, Fig. ^ also follow a remarkable pat- 
tern that differs from the statistical-model calculation. The latter predicts a 
monotonic, slowly rising trend with nuclear charge, due to Coulomb effects. 

It has been argued |]l8i that such obser- 




vations give strong evidence for a non- 
equilibrium situation, that can only be 
treated correctly by transport model 
approaches, such as Quantum Molecu- 
lar Dynamics (QMD).[0 It was also 
proposed [§] that a sizeable part of the 
kinetic energy stems from fast freeze- 
out of Fermionic motion. EOl . Whatever 
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Figure 22: Average kinetic energies of spec- 
tator fragments as function of their charge 
in the reaction Au+Au at lA GeV. SoHd 
symbols: data, open symbols: SMM calcu- 
lations. B 



the final explanation will be, such 'non- 
thermal' or 'non-classical' phenomena 
make the determination of the so-called 
caloric curve i.e. a plot of the 
temperature versus the excitation en- 
ergy per nucleon, very difficult. If 

^ is to rep- 



such a curve, see Fig. | 
resent a state property of (finite) nu- 
clear systems, perhaps suggesting by 
its plateau-like behaviour a liquid-to- 
gas transition, then it should be stable 
against the method applied. The major problem is to determine the abscissa 
which is usually obtained by the 'calorimetric' method, consisting essentially 
in adding up all the measured kinetic energies, perhaps with some high-energy 
tails cut off in a more or less well justified way. A preliminary result |p ob- 
tained recently. Fig. 0, shows that the apparent plateau width depends on 
the incident energy. 

Despite these unsolved problems in deriving a 'robust' caloric curve, it is 
highly interesting that theoretical approaches using fermionic dynamics [^, 
predict caloric curves that have great similarity to the curve presented in 

Fig. m 
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Figure 23: The caloric curve showing the 
dependence of the isotopic temperature on 
the excitation energy per nucleon. From 
ref. § 



Figure 24: Comparison of the caloric curves 
obtained for Au+Au at 600A and lOOOA 
MeV. From ref. S 
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Central collisions: 
Flow and system-size dependences 

The fireball is not enclosed in a box 
(thermal model) but is expanding in 
much the same way as the universe. At 
freeze out the configuration of nucleons 
is expected to look roughly as shown in 
Fig. There is a Woods-Saxon like 
density profile with a bulk part, having 
subsaturation density and an extended 
even lower-density surface part. Within 
the bulk part the flow velocity is ris- 
ing linearly. With such a scenario, and 
assuming that the local temperatures 
are not varying much over the configu- 
ration, one can describe kinetic energy 
spectra of emitted clusterized nucle- 
ons (Fig. |26|, left). A characteristic feature is that in contrast to a pure Boltz- 
mann scenario (no fiow) one sees curved mass dependent spectra. Those of 
heavier particles show a 'blueshift' in a spectacular way as they are less sen- 
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Figure 25: Freeze out profile (schematic) 
density and flow velocity (straight line) 
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Figure 26: Left panel: Kinetic energy spectra of Z—1 and Z=4 fragments in central Au on 
Au collisions. Right panel: Mass dependence of average kinetic energies. The smooth curves 
in both panels are blast-model fits to the data. 



sitive to thermal fluctuations that sit on top of the flow patterns and tend to 
wash out the flow features. Another way of showing evidence for radial flow are 
the mass-dependent average kinetic energies (Fig. |2^, right). In a thermalized 
system enclosed in a box, rather than expanding, the kinetic energies should 
be independent on mass. 

The energy taken up by such an expansive flow was estimated to be as high 
as 60%. 3 At 250A MeV incident energy, to be specific, this still meant that 
about 20 MeV per nucleon (the Q-values are about — 5A MeV) were available 
for 'local use'. It is therefore not surprising that statistical multifragmentation 
calculations, allowing for the available local energy, still failed to reproduce the 
measured fragment yields as illustrated in Fig. To obtain fair agreement 
with the data one has to assume that the primary fragments are cold ^3[, in 
contradiction to statistical expectations. 

Transition to transport models, not limited by the local equilibrium assump- 
tion, seems mandatory. In the past the IQMD model calculations [|T^ 
have been unable, though, to reproduce the observed IMF multiplicity, see 
Fig. E^. Modern developments are promising: a new cluster finding algo- 



rithm p5[ seems to be successful. 

One has tried to overcome the lack of quantum fluctuations in these quasi- 
classical simulations by introducing either so-called Quantum-Langevin meth- 



ods or by using Antisymmetrized Molecular Dynamics ^8 



In the following we discuss new, partially preliminary results obtained by 
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the FOPI Collaboration, that shed additional light onto this complex subject. 

It is of high interest to study system size dependences. Are the fireball 
data as 'universal' as the 'spectator' data are? We have studied the number of 
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tal (FOPI) charge distributions with hy- 
drothermal estimate 
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sions with reduced impact parameter below 
0.3. 



Figure 30: Same as adjacent figure, but for 
pions 
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emitted intermediate-mass fragments and the pion multiplicities varying the 
system size from ^^Ni+^^Ni, via ^^*^Xe+CsI to Au+Au, see Figs. ^ and 
Roughly, normalizing to a total of hundred observed baryonic charges, we find 
that three IMF are emitted in Au on Au reactions, two IMF in Xe+Csl and one 
IMF in Ni on Ni reactions. The order is inversed for pion emission: relatively 
more pions are emitted from the lighter Ni+Ni system (Fig. |5D|) . Using thermal 
language, one might say that heavier systems appear to be cooler at freeze-out, 
perhaps favoured by a more substantial bulk expansion. 

This seems to be corroborated by flow studies. Coming back to the direc- 
tivity D introduced in section 2 to sort out very central collisions, we compare 
in Fig. 31 directivity distributions measured for Ni+Ni and Au+Au at 90 and 
400A MeV. The reduced impact parameter range is < < 0.4 (touching 
configurations corresponding to h^^^ = 1). The (two-dimensional) random walk 
distributions of D (crosses), obtained by randomizing the azimuthal emission 
angles in each event, are compared with the actual distributions (histograms). 




Figure 31: Directivity distributions in various indicated systems at the same normalized 
centrality. The signal (histogram) is compared to distributions from azimuthally randomized 
events. 

One finds that the random distributions of this scale-invariant observable 
are broader for lighter than heavier systems (i.e. the relative fluctuations are 
larger), see the top two panels in Fig. and narrow down at higher energies 
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400A MeV 



0.15 



because of increasing particle multiplicities (compare the top with the bottom 
panels). Comparing with the actual distributions one sees that 

1) At 90 A MeV there is no net effect 
for the Ni system, while a statistically 
significant signal is seen for Au. 

2) At 400A MeV sideflow has increased 
substantially in both systems. 

3) The sideflow, Fig. 32, after correction 
for the finite number fluctuation, shows 
a dramatic system-size dependence. 

We are currently also studying other 
aspects of flow, such as radial flow and 
out-of-plane flow as a function of sys- 
tem size. It is intriguing to connect in- 
formation on flow with the probability 
to form clusters. The failure of the hy- 
drothermal attempt to understand the 
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degree of clusterization. Fig. | 

gests that rapid expansion physics may 
well require the consideration of fast non-equilibrium mechanisms, as proposed 
elsewhere |^ in a different context. The very low 'temperatures' suggested by 
the analyses |]T2|, |^ , see Figs. [1^, [2^, would then be an artifact of an inadequate 
assumption (equilibrium) . 



8 Non-equilibrium: 

isospin tracer method 

Valuable insight into the question of equilibration can be gained by studying 
systems, that, while mass symmetric, are isospin asymmetric. We have per- 
formed an 'isopspin tracer experiment' using four mass 96+96 systems Ru-|-Ru, 
Zr-|-Zr, Ru+Zr and Zr+Ru at 400A MeV. Naively, if Ru was to impinge on Zr, 
then, after the collision and in the case of complete transparency, we should 
be able to count 44 protons (Ru) in the forward rapidity hemisphere, and we 
should find the 40 protons of Zr at backward rapidities. In the extreme limit 
of a collective 'rebound' (predicted by one-fluid non-viscous hydrodynamics, 
as we shall see), just the opposite would be observed, while 'complete mixing' 
would show up as a 'white' distribution. 



A more complete description of this experiment is given in reference. |]30 
In the studied reactions, if central collisions are selected, the transverse- 
momentum vs rapidity distributions of emitted protons indicate the apparent 
dominance of one nearly isotropic source located at midrapidity. Fig. 33 shows 
for very central collisions the result of 'proton-counting' (including those in 
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Figure 33: Proton counting as a function 
of polar angle in four mass-96 on niass-96 
systems at 400A MeV (see text). 



Figure 34: Normalized triton/'^He ratios as 
a function of centrality 



deuterons) for the four systems as a function of the center-of-mass polar angle. 
To quantify conveniently the 'degree of mixing' we first study the symmetric 
systems Ru+Ru and Zr+Zr and renormalize in each angular bin the observ- 
able to +1 (Zr), resp. —1 (Ru). Using this normalization, we then study the 
asymmetric systems Ru Zr and the 'inverse' reaction Zr Ru. As can 
be seen, the two asymmetric systems converge to the expected mixed value 
(zero) only near 90°, where they must converge. As soon as one leaves the 
90° interval, the proton counting observable tends to the value expected for 
transparency, the sign being opposite for the Ru+Zr and Zr+Ru systems, as 
it should be for internal consistency. 

An alternative way of demonstrating this partial transparency is to use 
instead the triton to ^He ratio (Fig. 34). Again one proceeds by looking at 
the symmetric system first, normalizes the difference to ±1 and then studies 
the asymmetric systems. This time we fix the phase-space region: we look 
at to y^^^ = 0.25 — 0.75. We now vary the centrality from impact parameters 
below 1 fm (ER6) to about 8 fm (M2). Note again that effects can only be 
expected away from midrapidity. We find that our second observable again 
indicates partial transparency, although the effect is least pronounced for the 
most central collisions. One can show by the way ||3^ that the 'fully mixed' 
value should be somewhat below the zero line in Fig. 34, near (-0.20). 

If one uses the proton-counting method as a function of rapidity, rather than 
polar angle as was done in Fig. 33, one can actually decompose the rapidity 
distribution in a projectile and a target 'remnant', as shown in Fig. 35. 

What can we conclude? We believe these types of data are important for 
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Figure 35: Decomposition of the rapidity Figure 36: Passing time, Fermi time and 
distribution in central collisions of ^^Ru + mixing time (symbols) as a function of 
9f^Zr at 400A MeV. After ref. |3^ beam energy 



at least two reasons. 

1) it is good to have a signature of partial memory of the past: this puts a 
constraint on dynamic models; these models must show that they can realisti- 
cally cope with the observed partial transparency, independently whether its 
origin is connected with the surface corona and possibly even with insufficient 
stopping volume. Reactions with heavy ions always concern finite sizes and 
corona effects. 

2) If the system is not fully equilibrated at the end of the evolution, we expect 
it to be even less equilibrated at earlier stages, i.e. at maximum compres- 
sion. Hence the pressure expected from the nuclear equation of state (EOS) 
is not fully developed, which in turn will lead to decreased fiow at the end 
of the evolution. Transparency may mock up a 'soft' EOS. Transport models 
used to extract information on the EOS must therefore correctly reproduce 
the observed partial transparency (be it bulk and/or surface transparency). 

Clearly, corona effects must be one of the reasons for observing partial 
transparency, although they are subtle: for the most central collisions they 
do not appear as an obvious structure close to target or projectile rapidity 
(see Fig. 35). Various transport calculations shedding light on the mixing 
problem have already been done. In Fig. 36 we show results deduced from 
recent work of Hombach et al. Using a periodic box condition the authors 



have determined the equilibration time for nucleon mixing in the SIS energy 
range. They find times on the order of 40 to 10 fm/c decreasing with increasing 
incident energy (see the Figure). 

Comparing this with the passing time {2R/u) {R projectile radius, u four- 
velocity), and the Fermi-time (typical time for a nucleon to cross the system 
due to its Fermi- motion) , we can conclude from this study that the passing 
time might be somewhat too short for equilibration below 500A MeV. The 
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Figure 37: Momentum space distributions in a transparency and a hydro-shock scenario. 



IQMD - b < 1 fm 



relative transparency at lower energies is intimately connected with the im- 
portance of Pauli blocking and the still frozen nucleonic degrees of freedom. 
(The mobility implied by the Fermi time is important for the communication 
between 'participants' and 'spectators') . 

Transparency is sensitive to the 
medium- modified nucleon-nucleon elas- 
tic scattering cross sections (T„„. Fig. 37 
shows momentum-space distributions 
for Ru+Ru at 400 AMeV obtained us- 
ing the IQMD transport code with 
two extreme, but instructive, assump- 
tions: for the left panel we used 
0.5 * cr„„ and for the right panel 5 * cr„„ 
{arm are free-scattering cross sections 
from the literature). We note that the 
code of course applies the usual Pauli- 
blocking weights in the collision term 
[ P^ . This takes care of some, but not all 
in-medium effects. The two (impact pa- 
rameter zero) cases illustrate the trans- 
parent and the 'hydro-shock' alterna- 
tives. In the former case one still sees 
'spectator' sources, in contrast to ex- 
periment. In the latter case, almost instant equilibration leads to highly 
oblate event shapes with increased particle emission at 90° (c.o.m.): the true 
'squeeze-out' predicted in the early seventies |Q. The consequences for our 
proton counting observable are seen in Fig. 38. In particular the hydro-shock 
scenario leads to a rebound, the opposite of what we observe. The data are 
in between the two scenarios (although closer to the 'transparent' scenario) 




Figure 38: Comparison of the proton count- 
ing observable with various IQMD simula- 
tions 
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and are in impressive agreement with IQMD predictions using the nominal 
(sohd curve). It is important to reahze that only a fraction of the flow- 
creating pressure expected from ideal hydrodynamics is achieved in a heavy 
ion collision. 



9 Conclusion 

We have presented multifragmentation data from heavy ion reactions at inci- 
dent energies varying from 0.1 A to lA GeV. Two powerful and to some degree 
complementary experimental setups, ALADIN and FOPI, allowed us to study 
the complex decay of both 'spectator' and 'flreball' or 'participant' matter. 
It appears that rather different mechanisms govern these multiparticle events. 
While a remarkable 'universality', i.e. invariance under variation of the sys- 
tem size and of the incident energy, was seen to characterize the decay of the 
intermediate mass fragment source in the case of spectator matter, we saw 
that in contrast, the decay of rapidly expanding fireball matter showed pro- 
nounced system-size dependence, that could be traced, at least in part, to the 
occurrence of significant size- dependent fiow. 

Purely statistical concepts, comparatively successful in explaining the yields 
of particles with Z > 2, fail when it comes to explaining light particle yields, 
and above all when momentum space distributions are to be understood. Even 
the 'universal' spectator decays are accompanied by an apparently fast removal 
of some of the available energy primarily by single nucleons and show a non- 
trivial evolution of average kinetic energies with the nuclear charge, implying 
the presence of either fiow (sideways bounce or radial) or fast freeze-out of 
Fermionic motion, possibly due to large scale fiuctuations caused by the system 
passing through the spinodal region, or some entrance-channel memory effects. 

Such phenomena must be understood before we can think of deriving a 
'robust' caloric curve for finite nuclei, that we eventually could extrapolate to 
infinite nuclear matter. 

It seems also urgent to join up the 'spectator' data with information from 
lower energy reactions, both central (one-source) and peripheral (quasipro- 
jectile). Evidence for a liquid-gas transition has been claimed recently for 
quasiprojectile decay, as well as spectator decay. 0, ^ 



The mechanisms of multifragmentation in an expanding fireball are not 
yet fully understood. The underestimation of the degree of clusterization in 
a hydodynamic scenario with local equilibration, the observed system-size de- 
pendences and the lack of complete isospin-mixing indicate that fast expansion 
physics maybe a non-equilibrium process. While the role of fiow seems quali- 
tatively evident, the quantitative connection to the Equation of State must be 
further explored. 

New isospin-tracing experiments have a potential to overcome some of the 
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questions relating to a quantitative assessment of non-equilibrium features. 

The fact that we observe the clusterization of finite pieces of fermionic 
matter in at least a partially non-equilibrated environment calls for application 
and the further development of quantum transport theories, that fulfill in the 
static limit the well known requirements of nuclear physics theory. 
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